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One of the simplest proposed experimental probes of a Majorana bound-state is a quantized 
(2e 2 /h) value of zero-bias tunneling conductance. When temperature is somewhat larger than the 
intrinsic width of the Majorana peak, conductance is no longer quantized, but a zero-bias peak 
can remain. Such a non-quantized zero-bias peak has been recently reported for semiconducting 
nanowires with proximity induced superconductivity. In this paper we analyze the relation of the 
zero-bias peak to the presence of Majorana end-states, by simulating the tunneling conductance 
for multi-band wires with realistic amounts of disorder. We show that this system generically 
exhibits a (non-quantized) zero-bias peak even when the wire is topologically trivial and does not 
possess Majorana end- states. We make comparisons to recent experiments, and discuss the necessary 
requirements for confirming the existence of a Majorana state. 



Recent proposals [1-7] to build topological supercon- 
ductors from conventional spin-orbit coupled systems 
have sparked an active experimental effort to realize Ma- 
jorana fermions and probe their predicted non-Abelian 
exchange statistics. Tunneling from a normal wire into a 
topological superconducting wire with a Majorana end- 
state yields a quantized G(0) = ^- conductance peak 
at zero-bias [8-10]. This quantized zero-bias peak (ZBP) 
constitutes one of the simplest and most direct experi- 
mental probes for a Majorana fermion, and is likely to 
be the first test conducted on any putative topological su- 
perconducting wire. The observation of quantized zero- 
bias conductance with G(0) = ^f- ZBP requires tem- 
perature, T to be sufficiently smaller than the intrinsic 
width, 7, of the Majorana peak, due to hybridization 
with the normal lead. For T comparable to or somewhat 
larger than 7 a ZBP may still occur, but its no longer 
quantized and can take any value less than ^-[8-10]. 

A recent set of experiments on InSb nanowires coated 
with a superconducting NbTiN layer report the obser- 
vation of non-quantized ZBP's when a magnetic field of 
sufficient strength is applied along the wire [11]. Similar 
results have since been reported by other groups[12, 13]. 
These experimental observations are qualitatively con- 
sistent with the existence of Majorana end-states, and 
constitute an important first step towards the realiza- 
tion of Majorana fermions in solid-state systems. Given 
the potential significance of these findings, it is impor- 
tant to build a more quantitative understanding of the 
experimental system. In particular, we would like to es- 
tablish whether observed non-quantized ZBP's definitely 
correspond to thermally broadened peaks from Majorana 
end-states, or whether they could be produced by some 
other mechanism? 

To this end, we have conducted numerical simula- 
tions of tunneling conductance for spin-orbit coupled 
wires with proximity- induced superconductivity. Our 
simulations use realistic energy scales appropriate for 
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FIG. 1. (Top) Schematic of tunneling geometry. (Lower Left) 
Dispersion of sub-bands in multi-band wire. Each sub-band 
is split by fioB due to the magnetic field. Majorana fermions 
appear only when an odd number of sub-bands is occupied. 
(Lower Right) Color plot of tunneling conductance, G, at fi- 
nite temperature as a function of applied field /noB and lead- 
wire voltage, V, for a multi-band wire with realistic amounts 
of disorder (see Fig. 2 for detailed parameters). A stable zero- 
bias peak appears despite the fact that there is no Majorana 
end-state. At lower temperature, the peak is revealed to come 
from a cluster of low-energy states (see Fig. 2). 

InSb wires, and consider the various experimentally rel- 
evant non-idealities including: multiple occupied sub- 
bands[14-17], modest amounts of disorder[15, 18-20], 
and non-zero temperature. Our study reveals impor- 
tant features absent in previous studies of clean- or very 
weakly disordered wires [2 1-23]. 

We find that, at non-zero temperature and in the pres- 
ence of multiple sub-bands and weak disorder, zero-bias 
peaks generically occur even when the wire is in the topo- 
logically trivial regime and does not have Majorana end- 
states. Furthermore, we find that the ZBP's persist even 
when disorder is sufficiently strong to destroy the topo- 
logical phase and fuse the Majorana fermions on each 
side of the wire [15, 18-20]. Such ZBP's are also found 
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outside the range in chemical potential where Major ana 
end states are expected in the clean limit, and are pro- 
duced by ordinary fermion states localized to the wire 
ends and clustered near the Fermi-energy. These states 
are in some sense, remnants of Majoranas, and appear 
and disappear with magnetic field in the same way as 
true Majorana end-states. Therefore, we argue, that the 
only way to definitely rule out a non-topological origin 
to the ZBP is to lower temperature below the thermally 
broadened regime and observe a truly quantized zero-bias 
conductance peak, well isolated from other background 
states. 

Model — We consider a three-dimensional rectangular 
wire of length L x along the x direction and cross sectional 
area L y x L z in the yz-plane. The continuum Hamiltonian 
for the spin-orbit coupled wire without proximity induced 
superconductivity is: 
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Here aR is the Rashba velocity, related to the spin or- 
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coupling to the magnetic field B taken throughout to 
point along the wire (in the x direction), and Ao is the 
proximity-induced pairing amplitude. 

To model this system, we approximate the contin- 
uum Hamiltonian by the following discrete tight binding 
Hamiltonian, defined on a N x x N y x N z site prism: 
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Here, we have included a random impurity potential 
Vimp( r )> which is chosen independently for each site, iden- 
tically distributed according to a Gaussian with variance 
V(r)V(r') = W 2 5 r ,r'j where (• • • ) indicates averaging 
with respect to disorder configuration. 

Table I relates the tight-binding parameters to the 
continuum model and gives estimated values for InSb 
nanowires. There is considerable uncertainty in the esti- 
mated spin-orbit strength, which was measured without 
the superconducting layer [24]. Since E so derives solely 
from the inversion symmetry breaking potential of the 
substrate- wire and superconductor- wire interfaces, the 
actual value could be rather different than 50/ieV, and 
one should consider the possibility that E so is much 
smaller (or larger). The mean- free path from disor- 
der is i = vft where r _1 ~ 27r j^ L a iV(/x) is the elas- 
tic scattering rate due to impurities. Here N(/i) is the 
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FIG. 2. Conductance traces as a function of /io^, correspond- 
ing to the color plot in Fig. 1. From bottom to top fi^B 
ranges from to 3Ao in steps of 0.2Ao (curves are offset for 
clarity). \i — —172 corresponding to cut A in Fig. 3, and 
to 6 occupied sub-bands (including spin). ZBPs appear for 
fioB « Ao just as for a Majorana end-state, despite having 
an even number of occupied sub-bands. Wire dimensions are 
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FIG. 3. (Left) Color plot of conductance for same parameters 
as Fig. 2, but with fixed fioB = 1.5Ao- [i is varied so that 
between 5-6 sub-bands are occupied, and is measured with 
respect to the center of the 3rd topological region located at 
/i3 = —175.2. Stable ZBPs occur for a wide-range of /i, but 
have no topological origin. Detailed /JoB dependence is shown 
for fi corresponding to the dashed line A (see Figs. 1 and 2) 
and line B (see Fig. 5). The topological region in the absence 
of disorder, indicated by "Top. Region", does not exhibit 
ZBPs for this disorder configuration. 

ID density of states at the chemical potential: N(fi) = 



y 



where the sum is over occupied sub- 
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bands labeled by n and having band-bottoms located at 
energy e n . Transport experiments estimate t « 3/im 
(again without a superconducting layer) [25]. Since these 
measurements were done at large source-drain bias, this 
value reflects a sort of average over the lowest 3-4 sub- 
bands, and should be taken as a rough guide. 

Since only the outer-boundary of the wire is in contact 
with the superconductor, there will in general be different 
proximity induced gaps for different sub-bands. These 
multi-band effects can be important for reproducing the 
observed data for InSb wires. There, coherence peaks are 
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Parameter 


Symbol 


TB Equivalent InSb Value 




Wire Diameter 




N y , z a lOOnm 




Wire Length 


L x 


N x a ~ 2-3/xm 




Band Mass 


m 


(2ta 2 )~ 1 0.015ra e 


? 


Spin-Orbit 


E so - \mo? R 


Ur = VK^i 50/xeV 




Induced SC Gap 


A 


A 250/xeV[24] 


o 


Mean Free Path 


i 


(see text) « 3/xm[25] 




Min. Temperature 


T 


0.03A 





TABLE I. Tight-binding (TB) model parameters, and esti- 
mated values for InSb/NbTiN experiment [11]. a denotes lat- 
tice spacing in the TB model. 



observed at energy, Aq ~ 250/ieV, but non-zero conduc- 
tance occurs within the proximity induced "gap". The 
shape of this sub-gap conductance is not consistent with 
a Lorenzian broadening of the coherence peaks from cou- 
pling to the lead, but could be explained with a multi-gap 
scenario in which some sub-bands have proximity induced 
gaps smaller than Aq. However, we have so far not been 
able to reproduce, in detail, the large sub-gap tunneling 
conductance at ji^B = observed in the experiment. 
Tunneling Conductance - We use the iterative trans- 
fer matrix method to construct the Green's function for 
the end of the wire, and compute the scattering matrix 
from the Green's function[10, 26]. Experimentally, the 
wires are terminated by a large gap superconductor and 
only Andreev reflection contributes to tunneling current: 
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where r e h(s)ij is the electron- hole part of the reflection 
matrix from channel i to channel j in the lead, and / is 
the Fermi distribution. Throughout, we take the lead to 
have position and energy independent density of states 
No — 1/(ttvf), and model the tunneling barrier by a 
weak hopping link with hopping strength £lw- The lead- 
wire coupling is characterized by 7lw = ^o|^lw| 2 - 

Proximity induced superconductivity is modeled by 
coupling the boundary of each cross-section in the yz- 
plane to an infinite superconductor, producing the self- 
energy [18, 27]: 



^SC (<U+,y,z) = Vedge 
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The projection P e dge(y,z) = 1 if (y,z) is on the outer 
edge of the wire (y = 0) and is zero otherwise. Here, 
uj+ = uj + if] where r] positive and infinitesimal, T3 is the 
z-Pauli matrix in the Nambu/Gorkov particle-hole basis: 

^(k) = (c k)t c kj ; -cLk^ c^t) ' wnere C M destroys 
an electron with momentum k and spin s G {t,|}- 
is the pairing gap of the adjacent bulk superconductor 
7sc = 7tNb(0) |T| 2 is the strength of coupling between the 
wire and superconductor, Nb(0) is the density of states 




FIG. 4. Signature of an ideal Majorana state for a long, clean 
wire (N x = 360, W = 0). (Top Row) fioB dependence for \i 
in the center of the 6th topological region (11 occupied sub- 
bands), at T = (left) and T = 0.03 A (right). (Bottom) 
Corresponding color plot for T = 0. 



of the bulk superconductor near the Fermi-surface, and 
T is the wire-superconductor tunneling amplitude. 
Conditions for Majorana End-states — In a multi- 
band wire, when ji^B > Ao, one can think of each sub- 
band as contributing a Majorana end- state which then 
mix. For an even number of occupied sub-bands, the 
Major anas fuse into ordinary fermions and are pushed 
away from zero-energy. By contrast, for an odd num- 
ber of occupied sub-bands a single Majorana state al- 
ways remains at zero-energy [15]. More quantitatively, 
to observe a Majorana, \i must fall within specific in- 
tervals of size « ±\noB — Ao|. In [11] ZBP's are ob- 
served for hqB > Ao ~ 250/ieV. Indicating that Majo- 
rana end-states can exist only in narrow regions of /i of 
size « 250/ieV, or « ^ of the typical sub-band spacing 
E sh « 2.5meV. 

Problems can arise for short wires [21, 28]. If the 
wire is not sufficiently long, then the Majorana states 
on each end of the wire can overlap, splitting into or- 
dinary fermion states. For example, the InSb wires in 
[11], had aspect ratios of L x « 20 — 301^. Tight-binding 
simulations of such a wire, with the parameters as in 
Table I, show that Majorana end-states hybridize with 
splitting of a few percent of A even for clean wires. 
This issue is exacerbated by disorder, which is strongly 
pair-breaking [18-20] and reduces the gap protecting Ma- 
jorana end- states, allowing them to spread out and hy- 
bridize more strongly. Indeed, realistic amounts of disor- 
der corresponding to the estimated £ readily destroy the 
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FIG. 5. (Left) Angle dependence of non-topological ZBP is 
similar to that of a true Majorana derived ZBP (curves are 
offset for clarity). measures angle of B and the x-axis in 
the x-y plane. (Right) Color plot of /JoB dependence for \i 
corresponding to cut B in Fig. 3. 



quantized conductance peak from Majorana end-states. 

Despite these issues, ZBP's are still observed in InSb 
wires. This suggests that ZBP might occur rather gen- 
erally and may not necessarily tied to the presence of 
zero-energy Majorana end-states. Below, we will show 
that this is indeed the case, and that stable ZBP's fre- 
quently appear without Majorana states. 
Tunneling Data- To start, we first consider an idealized 
case, showing quantized Majorana peaks. Fig. 4 shows 
conductance data for a long, perfectly clean wire of length 
N x = 360 in the center of the 6th topological region (i.e. 
with 11 occupied sub-bands including spin). The left 
panel shows traces at T = 0. As /jlqB is increased the 
bulk pairing gap collapses, leaving behind a single quan- 
tized Majorana peak. At finite temperature, T = 0.03Ao, 
the narrow Majorana peak is thermally broadened and 
greatly reduced in height. Fig. 4 along with a color plot 
of the T = curves. We note that, even for a clean wire, 
the observation of this quantized peak requires much 2-3 
times longer wires than those used in [11], and an order 
of magnitude larger than in [13]. 

Under realistic experimental conditions for semicon- 
ducting wires, disorder is expected to play an impor- 
tant role[18-20]. The pair-breaking effects of disorder in 
a magnetic-field tend to collapse the proximity induced 
pairing gap, causing Majorana peaks to split away from 
zero-energy. Furthermore, disorder reduces the mini-gap 
to end-states from other occupied sub-bands. In the ab- 
sence of disorder, these end-states are pushed up above 
induced gap where they get absorbed by a continuum of 
extended states. In a field, disorder reduces the mini-gap 
splitting, causing these end-states to cluster near zero- 
energy. Our simulations show that, that clusters emerg- 
ing from mini-gap states frequently remain localized near 
the end of the wire, and at temperature larger than the 
mini-gap splitting merge into a single ZBP. 

Fig. 3 shows a color plot of tunneling conductance 
for fixed n$B = 1.5Aq and T = 0.03Aq, as a func- 
tion of \i near the 3rd topological region (with 5-6 sub- 
bands occupied). The modest amount of disorder in- 



cluded here (W = 12) corresponds to a very long mean- 
free path (JL « 10/im), but is nevertheless sufficient to 
destroy Majorana end-states. Despite this, zero-bias 
peaks frequently emerge (we note that these peaks ap- 
pear even more readily for stronger disorder, for example 
with W corresponding to £ ~ 2.5/im relevant to InSb 
wires). These peaks are typically stable over intervals of 
A/i « Ao although they are no-longer tied to the topo- 
logical region |/i — ^3 1 < yj (/hqB) 2 — Aq. Fig. 2 shows a 
typical example of this non-topological ZBP developing 
as a function of applied field, corresponding to \i along 
cut A in Fig. 3. At T = (left), one can resolve the 
ZBP into multiple peaks near zero-energy, however at 
T = 0.03Ao (right) the peaks are smeared into a single 
ZBP. 

These peaks appear and disappear under very similar 
magnetic field conditions as true Majorana states would: 
1) they appear only when /jlqB > A when pair-breaking 
effects or disorder become important, 2) they disappear 
when the field is rotated to point along the wire (see 
Fig. 5). However, unlike ZBP's tied to Majorana states, 
these non-topological ZBP's appear commonly through- 
out the range of chemical potentials between adjacent 
sets of sub-bands, meaning that they would be more read- 
ily observed without strongly tuning fi. To illustrate this 
last point, Fig. 5 shows a color plot of the B-field de- 
pendence of conductance for \i corresponding to slice B 
in Fig. 3. The results are similar to slice A, and further 
simulations show that this B-field dependence is generic 
whenever a ZBP is stable over an interval of \i. 

The features described here depend strongly on the 
size of the spin-orbit coupling. For weaker spin-orbit 
coupling, non-topological ZBPs appear even more read- 
ily and at even weaker levels of disorder. Conversely, if 
E so is substantially larger than estimated in [24], it be- 
comes difficult to reconcile the observations with a ZBP 
of non-topological origin. 

Discussion — We have shown that, under experimen- 
tally realistic conditions for semiconducting wires with 
rather modest amounts of disorder, Majorana end-states 
are destroyed and do not give rise to quantized ZBPs. 
Nevertheless, at finite temperature, ZBPs of a non- 
topological origin often appear due to the pair-breaking 
effects of disorder, which lead to clusters of low-energy 
states localized near the wire end. These non-topological 
ZBPs are typically stable with respect to variations of 
chemical potential and magnetic field, and appear and 
disappear under nearly identical conditions to those of 
true Majorana peaks. 

These results strike a note of caution for interpreting 
recent experimental evidence of Majorana states in tun- 
neling data[ll-13]. In order to truly identify Majorana 
end-states and rule out a non-topological origin of ob- 
served ZBPs, it is crucial to push to lower temperatures, 
and observe a quantized conductance peak. At suffi- 
ciently low temperature non-topological ZBPs will reveal 
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themselves as clusters of states, and can be distinguished. 
Spurious disorder induced peaks can also be ruled out by 
more complicated measurements, such as the 4tt Joseph- 
son effect [28], which does not survive once Majoranas are 
destroyed by disorder [29]. 

Our simulations indicate that to realize Majorana 
states in semiconductors, one likely needs to produce 
substantially longer and cleaner wires. These difficulties 
suggest that it may be beneficial to seek materials with 
larger spin-orbit [30]. 
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During the completion of this manuscript a related 
work appeared by Bagrets and Altland, which finds a 
disorder induced zero-bias peak in the disorder- averaged 
density of states that may be related to the ZBP dis- 
cussed here. However, we caution that since the tun- 
neling conductance into one point of the wire is not a 
self- averaging property, it may be problematic to make 
comparisons from disorder averaged properties. 
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